Measuring the mass of a sterile neutrino with a very short basehne reactor experiment 
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An analysis of the world's neutrino oscillation data, including sterile neutrinos, (M. Sorel, C. M. 
Conrad, and M. H. Shaevitz, Phys. Rev. D 70, 073004) found a peak in the allowed region at a 
mass-squared difference Am^ =0.9 eV^. We trace its origin to harmonic oscillations in the electron 
survival probability Vee as a function of L/E, the ratio of baseline to neutrino energy, as measured 
in the near detector of the Bugey experiment. We find a second occurrence for Am^ '^1.9 eV^. 
We point out that the phenomenon of harmonic oscillations of Vee as a function of L/E, as seen in 
the Bugey experiment, can be used to measure the mass-squared difference associated with a sterile 
neutrino in the range from a fraction of an eV^ to several eV^ (compatible with that indicated 
by the LSND experiment), as well as measure the amount of electron-sterile neutrino mixing. We 
observe that the experiment is independent, to lowest order, of the size of the reactor and suggest 
the possibility of a small reactor with a detector sitting at a very short baseline. 

PACS numbers: 14.60.Pq 

Keywords: neutrino oscillations, sterile neutrino 



Present data demonstrate that neutrinos change their 
flavor while propagating in vacuum and through matter. 
The evidence comes from solar neutrino experiments [l| , 
a long baseline reactor experiment atmospheric ex- 
periments 0] , and a long baseline accelerator experiment 
These experiments, together with the constraint im- 
posed by the CHOOZ reactor experiment Q, provide a 
quantitative Q determination of the mixing parameters 
and mass-squared differences for three neutrino oscilla- 
tions. Moreover, as the data becomes ever more precise, 
alternative explanations of the data are continually being 
ruled out 7,]. 

The lone datum that does not fit into the scenario of 
three neutrino mixing is the appearance result from the 
LSND experiment Q. An oscillation explanation of this 
result requires a neutrino mass-squared difference of at 
least 10~^ eV^ while the world's remaining data is com- 
patible with two mass-squared differences of the order of 
8 X 10-5 eV^ and 2 x 10"^ gy^. The addition of a sterile 
neutrino or neutrinos [9] has been proposed in an attempt 
to incorporate LSND into an analysis that would be con- 
sistent with the world's data. Other physical mechanisms 
have also been proposed as possible explanations 

Restricting our discussion to sterile neutrinos (with CP 
and CPT conserved), the simplest extension is the inclu- 
sion of a single sterile neutrino ll| . Such models fall into 
one of two classes. The 3-1-1 scheme adds the fourth neu- 
trino whose mass separation is much larger than the other 
three. In the 2-1-2 scheme, the LSND mass-squared dif- 
ference separates two pairs of neutrinos with the smaller 
mass-squared differences. Current analyses indicate that 
neither provides a compelling explanation of the data 
13, Hi- As four neutrinos do not seem to be sufficient, 
five neutrino scenarios have been investigated 13|, [14 1 . 

Here we assume the existence one sterile neutrino in 
a 3-1-1 scheme and propose a new way to experimentally 




FIG. 1: The value of Ax^ = - xlrin versus Am\ The 
solid (green) straight line is the result for a three-neutrino fit 
to data from Refs. P, S, i, i, i, i, [T!]. The solid (blue) 
curve is the result of a four-neutrino fit to this data, and the 
dashed (red) curve is the result if, in addition, we add Bugey 
[l5l |. The zero for the vertical scale (Xmin) is arbitrary. 



measure the associated mass-squared difference and mix- 
ing angle. For simplicity, we introduce only one sterile 
neutrino; however, our methods can be easily generalized 
to include additional sterile neutrinos. For the baselines 
and energies of interest here, only the mass-squared dif- 
ferences involving the extra neutrino will contribute to 
the oscillations. The mass of the three usual neutrinos 
can be taken as degenerate so that oscillations will be 
governed by only one mass-squared difference Am^. 

In Fig. 4 of Ref. there is a narrow peak in the 
allowed region which occurs at Am^ = 0.9 eV^. We can 
trace this peak to the Bugey reactor experiment [3], 
which is an electron anti-neutrino disappearance exper- 
iment with detector baselines of 15, 40, and 95 m. We 
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construct a model of the aforementioned neutrino exper- 
iments. This analysis 17[ produces mixing angles for 
three neutrino oscillations that are very similar to those 
of Ref. [1]. In Fig. [U we present versus the mass- 
squared difference, Am^, for three cases. 

The horizontal solid (green) line is the result of a three- 
neutrino analysis of the data from Refs. [E Sj S fl @1j 
LSND Q, and KARMEN [13]. The contains the no 
oscillation contribution from LSND and KARMEN. By 
definition the three neutrino results do not depend on 
Am^ yielding a straight line. The solid (blue) curve is 
the result of a 3-1-1 analysis. As expected, the nonzero 
LSND data utilizes Am^ and the additional mixing an- 
gles 014 and ^24 (the results are essentially independent 
of 0^4) to lower the x^. For Am^ < 0.03 eV^ the fourth 
neutrino does not contribute to the LSND or KARMEN 
experiment, and the reverts to the three neutrino re- 
sult as it must. The exact values of these two curves are 
irrelevant for the discussion at hand. We include them to 
provide a reasonable background for the dashed curve in 
which we add the Bugey experiment [l5| to the previous 
analysis, following exactly the analysis in Ref. 15]. For 
Am^ < 0.2 eV^, the dashed curve merges with the solid 
curve as Bugey does not contribute in this region. Be- 
yond this, the dashed curve contains fluctuations. This 
is because the parameters are in a region that runs along 
the edge of the Bugey excluded region where the is 
not smooth. This phenomena is also present in Ref. [l3| . 
These curves do not indicate the existence of a sterile 
neutrino. The addition of the CDHS [H, CCFR84J3], 
and NOMAD [23| experiments were found in Ref. [ijto 
largely offset the LSND indication of a sterile neutrino. 

The first important feature of the dashed curve is the 
narrow dip near 0.9 eV^. This corresponds to the peak in 
the allowed region found in Fig. 4 of Ref. [l3] . Note that 
we also find a second narrow dip near 1.9 eV^. This larger 
mass-squared difference is also present in Ref. [l3| : how- 
ever, in their analysis, the inclusion of other null result 
short baseline experiments jH, [3, [l^l almost completely 
suppresses this dip's significance. Clearly the narrow dip 
in the h ere and the narrow peak in the allowed region 
in Ref. [13| originates from the Bugey experiment. 

To understand the source of this dip, we turn to the 
in vacuo neutrino oscillation probability. Using a stan- 
dard extension of the MNS mixing matrix, the electron 
neutrino survival probability at these short baselines is 
approximately 



Vee{L/E) 9i l-sin2(2 6ii4)sin2 



(1) 



with (j) = 1.27Am'^L/E where Ato^ is in eV^, the base- 
line L is in m, and the energy E is in MeV. We can use 
this to determine which data in the Bugey experiment 
yield the dips in x^- Fig. [2] contains a plot of Vee versus 
L/E from the Bugey near detector located at i = 15 
m. The solid (red) curve is for the best fit parameters 
with a mass-squared difference of Ato^ = 0.9 eV^ and the 




L/E (m/MeV) 

FIG. 2: The electron survival probability Vee versus L/E. 
The data are from the near detector, L = 15 m, of the Bugey 
[isj experiment. The solid (red) curve results from a ster- 
ile neutrino associated with Am^ = 0.9 eV'^ and the dashed 
(violet) curve results from Am^ — 1.9 eV^. 



dashed (violet) curve for Am^ = 1.9 eV^. Although the 
statistics are poor, one can see that there is an harmonic 
oscillation in the data at frequencies which generate the 
dips in the x^ ■ The dips correspond roughly to a change 
in of ten. Whether these dips are real or statistical 
fluctuations is not the question. The point is that with 
improved statistics the existence of a sterile neutrino with 
Am^ in the appropriate range would produce a measur- 
able narrow minimum in the x^- The Bugey experiment 
is an existence proof for the validity of such an experi- 
ment. 

For a reactor, the neutrino spectrum and technology 
set the range of detectable neutrino energies to lie be- 
tween approximately 1 to 5 MeV. At a baseline of 15 m, 
the ratio L/E ranges from 3 to 15 m/MeV. Oscillations 
have undergone one cycle when (f> = tt. For a Am^ of 
0.9 (1.9) eV^, the oscillation length is 2.7 (1.3) m/MeV 
which results in approximately 5 (10) cycles in the al- 
lowed range. 

The amplitude of the oscillations for the Bugey data is 
about 1.25%. Twice this (peak to trough) is somewhat 
less than the average error bar thus giving the low sta- 
tistical significance to these dips. If the statistical error 
bars were one fourth this (sixteen times the data) an os- 
cillation pattern of this same magnitude would be very 
significant. This could be achieved by running longer 
and/or building a larger detector and/or using a more 
powerful reactor. 

Another consideration is the energy resolution, or the 
minimum size of the bin in E. In order to cleanly de- 
fine the oscillation length in L/E, four data points per 
cycle are needed, or a resolution of 0.7 (0.33) m/MeV 
for Ain^ — 0.9 (1.9) eV^. This corresponds to an en- 
ergy resolution of 10% (5%). For Am^ = 0.9 eV^, this 
resolution is less stringent than that of the Bugey exper- 
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iment by about 25%. For Am^ = 1.9 eV^, the needed 
resolution is about double that in the Bugey experiment 
thus requiring a total 32 fold increase in counts relative 
to Bugey. 

We estimated these numbers for an average energy neu- 
trino. The result will approximately hold for smaller L/ E 
as can be seen in Fig. [2l Focusing on an L/E below 5 
m/MeV, one could combine two bins into one wider bin, 
reducing the error bars so that they become comparable 
to the other error bars, while retaining the requisite num- 
ber of data points per oscillation. However, for the larger 
values of L/E equal spacing in E yields wider spacing in 
L/E, and the ideal spacing is not achieved. 

There is an absolute maximum mass-squared difference 
that can be reached by this type of experiment. This is 
set by the physical size of the reactor. If one gets to a re- 
gion where the length of a single oscillation is comparable 
to the dimensions of the reactor core, then neutrinos from 
the back of the reactor are incoherent with the neutrinos 
from the front of the reactor. Using a scale factor for the 
reactor core of 3 m and an average energy of 3.5 MeV, we 
find the maximum achievable mass-squared difference to 
be roughly 3 eV^. This number depends on the shape of 
the reactor core, the location of the detector with respect 
to the orientation of the cylindrical core, and the power 
distribution within the core. A straightforward calcula- 
tion for a given experiment is necessary to get more than 
a crude estimate. Unfortunately, this number is smaller 
than the 10 eV^ that is the lower end of a presently al- 
lowed region found in Ref. [13]. To reach this maximum 
sensitivity, an energy resolution of 3% is required. 

To determine the lower limit on the measurable mass- 
squared difference, we require that the oscillation phase 
reach at least </> = 7r/4. From this phase, one could de- 
termine the oscillation parameters without knowledge of 
the absolute flux. This gives a lower limit of Am^ — 0.05 
eV^ for L = 15 m. 

The measurement of Am^ is insensitive to the absolute 
normalization of the data; the oscillation length derived 
from the harmonic oscillations in the data determines 
Am^ . The amplitude of these oscillations determines di- 
rectly sin^(2 6'i4) which is also independent of the abso- 
lute flux if several oscillation cycles are measured. How 
this works out in the data analysis can be seen in Fig. [2l 
For small values of L/E, the data are coherent and thus 
the peak is quite near Vee = 1- The uncertainty in the 
norm of the data will necessarily be sufficient to allow the 
data to be uniformly adjusted such that the fit curve will 
have the peak for small L/E also quite near Vee = 1, the 
required physical value. A great advantage of the pro- 
posed experiment is that the normalization of the data, 
usually the largest systematic error in a neutrino oscilla- 
tion experiment, is nearly irrelevant. 

The above utilized a baseline L = 15 m. Doubling 
the flux by moving to L = 11 m would be attractive. 
We repeat the analysis for this value of L. The max- 



imum mass-squared difference which can be probed re- 
mains around 3 eV^ . The minimum and maximum value 
of L/E shift to 2.2 and 10.3 m/MeV, respectively. The 
midpoint for L/E moves to 5.1 m/MeV, and the mea- 
surement would cover 4 (7) cycles for Am^ = 0.9 (1.9) 
eV^. To reach 3 eV^, a resolution of 4% in the energy 
would be required. By lowering the maximum value of 
L/E the minimum sensitivity is raised, here to 0.06 eV^. 
Going to a smaller baseline increases the flux by L^, but 
the size of the reactor, hence its power, would necessar- 
ily decrease by L^. This produces a smaller number of 
cycles such that the number of energy bins decreases by 
an additional factor of L. The overall result is that the 
experiment, to lowest approximation, is independent of 
the size of the reactor. This is true if you are searching 
for the existence of a sterile neutrino. Since the accu- 
racy of the measured mass-squared difference would be 
increased by observing additional cycles, the last factor of 
L would not apply if the goal were a fixed error on Am^ . 
The larger width in the energy binning also reduces the 
resolution needed. For L = 11 m, the required resolu- 
tion increases to 6 (13)% for Am^ = 0.9 (1.9) eV^ For 
a discovery experiment, a small research reactor with a 
small detector sitting very near the core is an interesting 
option to consider. The detector might even be wrapped 
around the core to increase the count rate. 

The argument does require that the technology for neu- 
trino detection scales nicely as the size of the detector. 
The technical issues for doing this experiment and those 
involved with using short baseline neutrino detectors for 
nonproliferation monitoring 21| are related. 

We have shown that a very short baseline reactor ex- 
periment can be used to measure the mass-squared dif- 
ference associated with a sterile neutrino should it lie in 
the range of less than a tenth of an eV^ up to several eV^ 
by measuring the oscillations in the data over a number 
of oscillation lengths. A Fourier transform analysis of 
such data would be an efficient way of extracting the os- 
cillation frequencies; however, such methods are not fun- 
damentally different from fitting oscillation parameters 
to the data as a function of the mass-squared difference 
Am^. The experiment also measures sin^(2 0i4) through 
the amplitude of the oscillations. We suggest that a small 
reactor with a very short baseline be investigated. The 
experiment is not sensitive to the absolute normalization 
of the data and thus holds the possibility of being more 
accurate than alternatives. 

The MiniBooNE experiment will soon confirm or 
contradict the LSND experiment. The experiment pro- 
posed here could play a significant role independently 
of that outcome. Should MiniBooNE leave the situation 
ambiguous, the proposed experiment could provide a cost 
effective, accurate, and independent way to resolve the 
situation. Should MiniBooNE confirm the existence of 
a sterile neutrino (s), then this experiment might provide 
a more accurate measurement of the mass-squared dif- 
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ference and of 614. Should MiniBooNE not see evidence 
for a sterile neutrino, it would be setting upper limits on 
the mixing angles. This proposed experiment might then 
be a way of further looking for a sterile neutrino in this 
range or, should a null result be found, further reducing 
the allowed value of 614. 

This work is supported in part by a grant from the US 
Department of Energy DE-FG02-96ER40975. 
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